Quick iterative scheme for the calculation of transfer matrices : application to MO ( 100 )

The transfer matrix of a solid described by the stacking of principal layers is obtained by an iterative procedure which takes into account 2" layers after n iterations, in contrast to usual schemes where each iteration includes just one more layer. The Green function and density of states at the surface of the corresponding semi-infinite crystal are then given by well known formulae in terms of the transfer matrix. This method, especially convenient near singularities, is applied to the calculation of the spectral as well as the total densities of states for the (100) face of molybdenum. The Slater-Koster algorithm for the calculation of tight-binding parameters is used with a basis of nine orbitals per atom (4d, 5s, 5p). Surface states and resonances are first identified and then analysed into orbital components to find their dominant symmetry. Their evolution along the main symmetry lines of the two-dimensional Brillouin zone is given explicitly. The surfacestate peak just below the Fermi level (Swanson hump) is not obtained. This is traced to the difficulty in placing an appropriate boundary condition at the surface with the tight-binding parametrisation scheme.